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1.  INTRODUCTION 

The  interaction  between  inclusions  and  cracks  has  been  an  active  field  of  study  for  many  years. 
Applications  ot  the  results  of  theoretical  analyses  have  wide  implications  in  the  implementation  of  a 
broad  range  of  new  and  emerging  material  systems.  The  understanding  gained  by  the  solution  of 
appropriate  elasticity  problems  allows  researchers  to  gain  insight  into  the  mechanisms  of  strengthening 
and  toughening,  as  well  as  material  damage,  because  of  the  presence  of  material  defects.  The 
mechanisms  of  crack  growth  in  composite  materials  and  ceramics,  and  that  of  strain  hardening  in  metal 
alloys,  are  primary  examples  of  the  direct  application  of  the  results  of  the  study  of  crack  interaction. 

The  problem  of  the  interaction  between  a  circular  inclusion  and  a  crack  has  been  solved  by 
Atkinson1  for  a  radial  crack  and  by  Erodgan,  Gupta,  and  Ruwani2  for  an  arbitrarily  oriented  crack. 
Erdogan  and  Gupta3  later  solved  the  problem  in  which  the  crack  crosses  the  interface.  These  solutions 
are  based  on  the  Green’s  function  solutions  consisting  of  a  dislocation  interacting  with  a  circular 
inclusion  (Dundurs  and  Mura4  and  Dundurs  and  Sendeckyj).5  In  an  earlier  paper,  Santare  and  Keer,6 
presented  the  two-dimensional  solution  for  the  dislocation  outside  a  rigid  elliptical  inclusion.  The 
problem  was  solved  using  the  complex  potential  methods  of  Mushkhelishvili.7  Special  attention  was 
paid  to  the  rotation  of  the  elliptical  inclusion  and  the  effect  that  has  on  the  stress  field  around  the 
ellipse.  Results  were  compared  with  a  power  series  solution  found  by  Stagni  and  Lizzio,1  which  did 
not  take  into  account  the  rotation  of  the  ellipse.  The  comparison  showed  that  the  rotation  has  a 
significant  effect  on  the  stress  field  in  many  instances.  In  another  paper,  Santare,  Keer,  and  Lewis9 
solved  a  related  problem  of  an  elliptical  hole  at  a  bone/implant  interface,  with  symmetrical  cracks 
radiating  from  the  edges  of  the  ellipse,  along  the  x-axis.  Solutions  to  the  resulting  singular  integral 
equations  were  found  using  a  numerical  scheme  proposed  by  Gerasoulis.10 

This  paper  provides  the  numerical  solution  of  the  general  problem  of  a  straight  crack  near  a  rigid 
elliptical  inclusioa  The  problem  is  formulated  using  the  results  of  Reference  6  and  the  numerical 
technique  used  in  Reference  9.  Tire  ellipse  can  be  made  to  have  any  shape,  from  a  circular  inclusion 
to  a  line  inclusion.  The  crack  is  oriented  at  an  arbitrary  angle  to  the  x-axis,  with  the  length  of  the 
crack,  L  =  Izj  -  z,l.  The  geometry  of  the  problem  is  shown  in  Figure  1. 
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2.  METHOD  OF  SOLUTION 

The  problem  is  posed  in  terms  of  the  complex  potentials  of  Mushkhelishvili,7  in  which  stresses 
and  displacements  can  be  written  in  terms  of  two  analytic  functions  as  follows: 


°xx  +  °yy  =  2  [  (z)  +  <t>'  (z)  1 

(1) 

°yy  '  <*xx  +  2i  Oty  =  2  [  Z  (z)  +  \j/'  (z)  J 

(2) 

2p  (u+iv)  =  K  (|)  (z)  -  z  <t>'  (z)  -  \|t  (z) 

(3) 

where  z  =  x  +  iy,  4  is  the  shear  modulus  of  the  matrix  material,  k  =  3  +  4  v  for  plane  strain,  and  k 
=  (3  -  o)/(l  +  v)  for  plane  stress,  where  o  is  the  poisson’s  ratio  for  the  matrix.  The  prime  notation  is 
understood  as  the  derivative  with  respect  to  z,  and  the  superimposed  bar  denotes  the  complex 
conjugate. 


2 


The  geometry  of  the  problem  is  simplified  by  using  the  function  that  maps  the  region  outside  the 
ellipse  onto  the  region  outside  the  unit  circle,  as  shown  in  Figure  2.  This  mapping  can  be  written  as, 


a  +  b 


a  -  b 
a  +  b 

and  a  and  b  are  shown  in  Figure  1. 


where 


R  = 


m  = 


Z-Plane  W-Plane 
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The  solution  for  a  dislocation  interacting  with  a  rigid  elliptical  inclusion  is  used  as  a  Green’s 
function  for  the  crack  problem.  The  potentials  were  derived  in  Reference  6  and  are  restated  here  for 
convenience.  They  include  the  dislocation  in  the  form  of  y  =  p.  (bx  +  iby)/in  (k  +  1),  the  mapped 
coordinates  C  and  C0-  which  are  the  trar.sformed  z  and  z0,  and  the  rotation  of  the  ellipse  e0. 


<t>  (C)  =  —  tog  [R  (C-Co  +  m/C  -  m/Co)] 

K 

+  log  KC  -  1/CoVCl  -  Y  log  I(C  -  m/CoVCl 
Y 

von 


(1/Co  ~  m/C o)  (1/Co  -  Co) 
(1/Co  *  Co/m)  (1/Co  *  0 


+  2pi  e0  R  m/KC 


(5) 


V  (0  =  Y  log  IR(C  -  Co  +  m/C  -  m/Co)]  -  Y 


Co  +  m/Cp 
C  -  Co  +  m/C  -  m/Co 


+  K7  log  [(C  -  1/CoVCl  ‘  Y  log  KC  -  m/CoVCl 


(m/Cp  ~  1/Co)  (m/Co  ~  Co/m) 
(m/Co  *  Co)  (m/Co  -  C) 


+  2p  ieD  R/C 


-  c 


+ 


1  +  mC2  r J_  1/Co  y  m/Co 

C2-m  C(C-1/Co)  C  (C  *  m/Co) 

_L  (1/^°  '  Q  -  p£o  m/tcC2] 

““  0/Co  Co/m)  (C  -  1/Co)2 


(6) 
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where 


e0  =  Re  i£-  Ky/^  +  ym/C0  +  ym2/^ 

(m/Co  -  1/Co)  (mCo  '  Cc/m> 

-ym/K'Co  -  ym  - - - 

(m/Co  ‘  Co) 

+  _L  (1/Cq  '  ^Co)  UCo  -  Co)  I  /2(xr  n  +  m  )  (7> 

K  (Vlo-lo/n)  J  K 

The  first  term  in  <J)  and  the  first  wo  terms  in  y  represent  the  potentials  for  dislocation  in  an 
unbounded  medium.  These  terms,  when  applied  to  equations  (1)  and  (2)  to  calculate  the  stresses, 
become  singular,  as  l/\z-z0.  The  remainder  of  the  terms  in  both  potentials  account  for  the  interaction 
between  the  ellipse  and  the  dislocation.  These  terms,  when  used  in  equations  (1)  and  (2),  constitute 
the  nonsingular  stresses  due  to  the  interaction  between  the  dislocation  and  the  inclusioa  Note  also 
that  the  rotation  of  the  ellipse  is  taken  into  account  through  these  nonsingular  terms  and  is  called  e0. 
The  last  term  in  the  expression  for  y  (the  bracketed  terms)  contains  the  first  derivative  of  the 
interaction  portion  of  the  first  potential  <j).  This  is  a  natural  consequence  of  the  method  of 
Mushkhelishvili.7 

Summing  these  stresses  along  the  crack  length,  and  setting  them  equal  to  the  stresses  due  to  the 
external  load,  two  singular  integral  equations  result.  The  solution  of  these  equations  will  give  us  the 
distribution  of  dislocations  along  the  crack.  If  we  resolve  the  stresses  into  components  normal  to  the 
crack  (Mode  I),  and  along  the  crack  (Mode  II),  the  two  equations  can  be  stated  as  follows: 

(z,  zj  bn  (zj  d  z0  =  fn  (z)  (8) 


(z,  z0)  b,  (Zo)  d  z0  =  f,  (z)  (9) 


where  the  n  and  t  subscripts  refer  to  the  normal  and  tangential  components  with  respect  to  the  crack, 
and  fn  and  ft  are  the  stresses  due  to  the  external  load.  The  first  term  in  each  of  these  equations 
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represents  the  Cauchy  singular  portion  of  the  stresses  (as  noted  above),  and  the  second  term  contains 
the  nonsingular  parts.  The  normal  and  *  ..igential  cartesian  stress  components  can  be  determined  using 
a  standard  trigonometric  transformation.  Furthermore,  if  equations  (8)  and  (9)  are  rewritten  in  terms  of 
lz-z0l  e  one  can  determine  the  stress  due  to  the  Cauchy  singular  portion  of  the  above  integral 
equation,  in  terms  of  the  angle  0  (the  angle  of  the  crack),  as  follows: 


1 


°yy  = 

1  z-z0  1 

lb> 

(cos  0  +  cos  3  0)  -  bx 

(3  sin  0  +  sin  3  0)] 

GO) 

1 

ST  — 

(by 

(3  cos  0  -  cos  0)  -  bx 

(sin  3  0  -  sin  0)] 

(ID 

'-'xx  - 

l  z-z0  1 

1 

ST  - 

[by 

(sin  3  0  -  sin  0)  +  bx 

(cos  3  0  +  cos  0)] 

(12) 

wxy  - 

1  z-z0  1 

The  technique  used  to  solve  for  the  unknown  distribution  of  dislocations  was  developed  by 
Gerasoulis.10  He  uses  a  piecewise,  quadratic,  pulynominal  representation  of  the  singular  and 
nonsingular  parts  of  the  integral  equation  to  discretize  the  equation  into  a  set  of  algebraic  expressions 
suited  for  a  matrix-type  solution.  In  order  to  solve  the  equation,  we  must  first  make  an  assumption 
about  the  nature  of  the  singularity  at  the  tips  of  the  crack.  Because  the  crack  is  contained  entirely  in 
the  matrix  material,  the  stresses  at  both  crack  tips  will  be  singular,  as  l/Vz-z0.  The  dislocation 
density,  therefore,  will  also  exhibit  square-root  singularity  at  both  tips  and  can  be  rewritten  in  the 
following  form: 


b(t)  =  g(t)  /  VTri"  (13) 

where  t  +  2(z-z,)/(z2-z1)  -1.  The  functions  g(t)  are  continuous  and  bounded  on  the  interval  t  e 
[1,  -1).  For  a  given  n,  this  integral  is  broken  into  2n+l  integration  points  and  2n  collocation  points. 
Gerasoulis10  states  that  the  method  works  equally  well  for  unequal  meshes,  but  in  this  problem  an 
equal  spacing  of  the  integration  and  collocation  points  is  used.  Convergence  in  this  problem  is 
relatively  good  using  this  method.  A  value  of  n  above  6  was  unncessary  for  crack  lengths  on  the 
order  of  the  size  of  the  inclusion. 

3.  RESULTS  AND  CONCLUSIONS 

Figures  3  through  8  are  given  in  an  attempt  to  compare  and  contrast  the  present  results  with  those 
found  by  Atkinson1  and  Erdogar.,  et  al.,2  for  the  circular  inclusion  and  with  Grief  and  Sanders11  for  a 
line  inclusion.  In  the  last  of  these  papers,  the  integral  equation  is  formulated  somewhat  differently  and 
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Normalized  Stress  Intel 


corresponds  only  to  the  interaction  terms  in  the  current  problem.  The  present  solution  lakes  into 
account  both  the  interaction  terms  and  the  stress  singularities  associated  with  the  crack,  in  the  same 
form  as  Atkinson1  and  as  Erdogan  and  Gupta.3  In  Figures  3  through  7,  the  crack  is  radially  oriented 
with  respect  to  the  origin,  that  is  8  is  equal  to  a  Figure  8  presents  results  for  a  crack  in  which  0 
and  a  are  not  equal. 

Figure  3  is  a  plot  of  the  normalized  mode  I  stress  intensity  factor  Kj/o  °°Vl~ versus  the  distance 
from  the  close  tip  of  the  crack  to  the  inclusion  for  several  elliptical  shapes  (m  varying  from  0  to  1), 

with  the  angle  a  either  equal  to  0  or  90°  (along  the  x-axis  or  the  y-axis).  Note  that  the  stress 

intensity  factor  is  lowest  when  approaching  the  crack  from  the  side,  (a  =  90°)  as  would  be  expected. 
An  interesting  result  is  shown  in  Figure  4,  in  which  normalized  stress  intensity  factor  is  plotted  against 
a  for  varying  elliptical  aspect  ratios.  The  only  parameters  varied  in  this  figure  are  a  and  the  aspect 
ratio  of  the  ellipse.  The  distance  from  the  close  tip  of  the  crack  to  the  ellipse  and  the  crack  length 
are  kept  constant.  The  close  tip  of  the  crack  is  0.1/R  from  the  ellipse,  and  the  normalized  crack 
length  (L/R)  is  taken  as  1.  Note  that,  for  relatively  flat  ellipses  (ra  >  .5),  the  location  of  the  close  tip 
of  the  crack  varies  rapidly  with  small  changes  in  angle  a.  In  effect,  once  the  crack  rotates  slightly 
from  the  x-axis,  it  is  shielded  by  the  flat  side  of  the  ellipse,  and  the  stress  intensity  factor  drops 
rapidly.  Figures  5  and  6  illustrate  this  in  more  detail.  In  these  figures,  stress  intensity  factor  is 
plotted  against  distance  from  the  ellipse  for  small  a.  Figure  5  shows  the  results  for  an  ellipse  with  an 

aspect  ratio  of  0.8,  and  Figure  6  shows  the  same  for  an  ellipse  of  0.9  aspect  ratio. 

In  Figure  7,  a  different  approach  is  taken.  In  this  figure,  the  close  tip  of  the  crack  is  a  fixed 
distance  from  the  origin,  rather  than  from  the  ellipse,  and  the  angle  a  is  varied  from  0  to  20°  for 
different  aspect  ratios  of  the  ellipse.  The  normalized  distance  from  the  origin  to  the  close  tip  of  the 
crack  in  this  figure  is  2.1/R.  For  an  a  of  0°,  and  an  ellipse  with  an  aspect  ratio  of  1,  this  is 
equivalent  to  having  the  crack  0.1/R  from  the  ellipse.  Note  in  Figure  7  that,  for  all  angles,  the  stress 
intensity  factor  is  very  close  to  unity.  This  is  contrasted  with  about  a  25  to  30  percent  decrease  of 
stress  intensity  factor  for  these  small  values  of  a  when  the  crack  is  shielded  by  the  ellipse,  as  shown 
in  Figures  5  and  6. 

This  effect  has  broad  implications  is  the  modeling  of  toughening  mechanisms  in  brittle  materials. 

It  seems  that  a  crack  at  the  end  of  a  small,  thin  inclusion  would  be  likely  to  grow,  whereas  a  crack 
oriented  normal  to  the  major  axis  of  an  elliptical  inclusion  would  be  less  likely  to  grow.  Orientation 
of  the  inclusions-as  well  as  their  shape-must,  therefore,  be  taken  into  account  in  the  toughening 
models. 
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Figure  8  presents  results  for  a  crack  that  is  not  radially  oriented.  In  this  figure,  the  crack  is 
oriented  either  horizontally  (8  =  0°)  or  vertically  (0  =  90°),  while  the  y-coordinate  of  the  close  tip  of 
the  crack  (zj)  is  kept  at  a  value  of  0.5.  The  crack  length  is  maintained  at  a  constant  value  of  1.0/R, 
and  the  x-coordinate  of  the  close  tip  of  the  crack  (z,)  is  varied  from  0  to  a  value  of  3.0/R,  for  an 
elliptical  aspect  ratio  of  0.8.  Note  that  the  crack  becomes  unshielded  at  a  distance  of  approximately 
1.80/R  from  the  y-axis.  With  the  crack  oriented  horizontally,  the  stress  intensity  factor  dips  slightly, 
as  the  close  tip  of  the  crack  (zj)  approaches  the  end  of  the  ellipse.  For  a  vertical  crack,  the  stress 
intensity  factor  passes  near  unity  at  this  point  and  rises  to  a  value  significantly  above  1.  The 
interaction  with  the  ellipse  is  still  evident  for  this  crack  at  a  large  distance  from  the  ellipse. 

The  technique  presented  in  this  paper  is  extremely  flexible  in  nature,  as  there  are  several 
parameters  that  can  be  varied  independently,  and  in  combination  with  each  other,  to  produce  interesting 
and  useful  results.  Consequently,  there  is  an  infinite  number  of  possibilities  that  could  be  explored. 
This  paper  has  presented  just  a  few  of  those  possible  results.  The  parameters  varied  in  this  study  were 
chosen  in  an  attempt  to  compare  this  formulation  with  previously  published  results. 
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